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ABSTRACT 

We study the mass distribution of a sample of 28 galaxy clusters using strong and 
weak lensing observations. The clusters are selected via their strong lensing prop- 
erties as part of the Sloan Giant Arcs Survey (SGAS) from the Sloan Digital Sky 
Survey (SDSS). Mass modelling of the strong lensing information from the giant arcs 
is combined with weak lensing measurements from deep Subaru/ Suprime-cam images 
to primarily obtain robust constraints on the concentration parameter and the shape 
of the mass distribution. We find that the concentration Cvir is a steep function of 
the mass, Cyir oc 

^-o.59±o.i2^ with the value roughly consistent with the lensing-bias- 
corrected theoretical expectation for high mass (~ IO^^/i^^Mq) clusters. However, the 
observationally inferred concentration parameters appear to be much higher at lower 
masses (~ lO^^/z^^M©), possibly a consequence of the modification to the inner den- 
sity profiles provided by baryon cooling. The steep mass-concentration relation is also 
supported from direct stacking analysis of the tangential shear profiles. In addition, 
we explore the two-dimensional shape of the projected mass distribution by stacking 
weak lensing shear maps of individual clusters with prior information on the position 
angle from strong lens modelling, and find significant evidence for a large mean el- 
lipticity with the best-fit value of (e) = 0.47 ± 0.06 for the mass distribution of the 
stacked sample. We find that the luminous cluster member galaxy distribution traces 
the overall mass distribution very well, although the distribution of fainter cluster 
galaxies appears to be more extended than the total mass. 

Key words: dark matter — galaxies: clusters: general — gravitational lensing 



* Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. Based on 
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Inc., under a cooperative agreement with the NSF on behalf of 
the Gemini partnership: the National Science Foundation (United 
States), the Science and Technology Facilities Gouncil (United 
Kingdom), the National Research Gouncil (Ganada), GONIGYT 



1 INTRODUCTION 

Gravitational lensing plays a dominant role in determin- 
ing the mass distribution of distant galaxies and clusters of 
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galaxies, because it allows direct measurements of the dis- 
tribution of dark matter which accounts for ^ 90% of the 
total matter content of the Universe. Precise measurements 
of dark matter distributions are important not only for un- 
derstanding the formation of galaxies in the context of the 
hierarchical structure formation scenario, but also for test- 
ing the properties of the putative dark matter particle, in 
particular its cold and collisionless nature. Observations of 
clusters are suitable and very apt for the latter purpose, as 
the long cooling timescale for hot gas in clusters indicates 
that the cluster gravitational potential is mainly determined 
by the dynamics of the dominant dark matter component 
which are well predicted by A^-body simulations. However, 
it is also expected that the baryonic component in clusters, 
which dominates the mass in the innermost regions, should 
play an important role at some point. 

There are two regimes of gravitational lensing that can 
be employed to measure and map the mass distribution in 
clusters. One is strong lensing, i.e., the regime that produces 
drastic lensing events with highly elongat ed arcs or mul- 
tiple images of background objects (e.g., IKochanekl l2006l : 
iKneib &: Natarai an 2011, for reviews). The other lensing 
regime is weak lensing wherein statistical measurements 
of small distortions in the shapes of background galax- 
ies are produced by a massiv e, foreground cluster (e.g., 
iBartelmann Schneiderl 12001. for a review). Since strong 
and weak lensing probe mass distributions at different radii, 
the combination of these two is powerful and essential for 
the full understanding and detailed mapping of the gravita- 
tional potential of clusters. 

Lensing studies of clusters have indeed confirmed sev- 
eral important predictions of the standard A-dominated 
Cold Dark Matter (ACDM) model. For instance, from nu- 
merical simulations, the radial run of the density of the 
mass distribution is predicted to have a uni versal form , 
that of the Navarro- Frenk- White profile fe.g. jNavarro et al.l 
Il996l . ll997l l. Detailed lensing measurements of massive clus- 
ters have confirmed this prediction by and large, although 
there does appear to be some cluster-to-clu ster variations 
(e.g..lOkabe et al.l l2010l : lUmetsu et al.l l2011b). Another im- 
portant prediction is that massive clusters are on average 
highly non- sp herical with a typ ical major to minor axis 
ratio of ^ 2:1 (jjing Sutoll2002l l. Weak lensing measure- 
ments of dark matter distributions in a sample of mas- 
sive clusters have directly confirmed this prediction as well 
(|Oguri et al.ll2010l ). Furthermore, lensing analysis of several 
merging clusters of galaxies strongly suppor ts the collision- 
less nature of dark matter (Clowe et al. 2006: 1 Jee et aLlfioO?! : 
'Bradac et al. 2008; Merten et al.ll201in . although a possible 
case against the collis ionless nature has also been discovered 
(jMahdavi et al.ll2007l ). 

However, measurements of the concentration parame- 
ter, which appears in the Navarro- Frenk- White profile, from 
gravitational lensing observations for clusters have been con- 
troversial. The concentration parameter is defined as the 
ratio of the virial radius to the scale radius, and is a di- 
mensionless parameter that quantifies the degree of the 
mass accumulation in the innermost regions. The standard 
ACDM model makes a prediction about this parameter too, 
namely that more massive haloes or haloes at higher red- 
shifts have s maller values of the co ncentration parameter on 
average fe.g.. iBullock et 'aDl200ll ). However, measurements 



of the dark matter distribution in the massive lensing cluster 
A 1689 for instance, using both strong and weak lensing, have 
indicated that the mass distribution is surprisingly highly 
concentrated with an estimated concentration parameter of 
Cvir ^ 12, which is significantly higher t han the expected 
standard ACDM prediction of Cy ir ^ 4 (jBroadhurst et al] 
I2OO5I : lUmetsu Broadhurstll2008l l. 

The interpretation of this result of lensing clusters 
being over- concentrated requires careful consideration due 
to the role played by projection and selection effects. 
This is because lensing observables can only measure 
the projected mass distribution, and hence the recovery 
of the three-dimensional mass distribution requires addi- 
tional assumptions about the elongation along the line- 
of-sight direction. For instance, because of the large tri- 
axiality of dark matter haloes, the mass and concen- 
tration parameter inferred from gravitational lensing de- 
pend strongly on the orientation with respect to the line- 
of-sight direction. This naturally implies that both the 
mass and concentration are signifi cantly overestimated when 
obser v ed along th e majo r axis ([Clowe. De Lucia, King 
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The large orientation dependence in turn sug- 
gests that strong lensing selected clusters represent a 
highly biased population with their major axes preferen- 
tially aligned with the line-of-sight directio n (iHennawi et al.l 
I2OO7I : .O guri Blandfo rd 2009; Mene ghetti e t al.l \2{)1^. 
Thus, the dark matter distribution of A1689 alone 
does not pose a severe challenge to the ACD M model 
(Oguri et al. 2005; Oguri^ Blandford 2009; Cor less et al] 
2009; Sereno, Jetzer, Lubinil 120101 : ICoe et all I2OIOI: 
Morandi. Pedersen. LiniousinI I2OIII : ISereno Umetsul 
boill : I Morandi et alj l201lh , although such high values of 
the concentration parameter appear to be common in the 
combined stro ng and weak lensing analysis of massive 
clust e rs (IComerford &: NataraianI l2007l: iBroadhurst et al 



20081: lOguri et al.1 l2009l : lUmetsu et al.1 l2011bl : IZitrin et al 



201 Ibh . High concentrations can be tested with the distri- 
bution of the Einstein radii, but the res ults have still not 
converged (Broadhurst & Barkana 2008; Oguri k, Blandford 
| 2009l: [Richard et al.1 12OIOI : ICralla et al.1 12OIII : IZitrin et alJ 
l2011al la). Also puzzling is the fact that such high concen- 
trations have not been claimed in we ak lensing analysis 
("Okabe et al.' 2010) or X-ray analysis (Buot e et all l2007l : 
Ettori et al. 2010) for samples of massive clusters. 

While the combined strong and weak lensing analysis al- 
lows accurate and robust measurements of the concentration 
parameter, the current main limitation is the small num- 
ber of clusters available for such a detailed combined an aly- 
sis. The Sloan Ciant Arc s Survey (SGAS; Hennawi et alJ 
I2OO8I : iBavhss et all l2011al ). which is a survey of strongly 
lensed giant arcs from the Sloan Digital Sky Survey (SDSS; 
York et al.i I2OOOI ). has already discovered more than 30 
bright giant arcs and therefore offers an ideal technique to 
expand the sample of clusters appropriate for detailed lens- 
ing analysis. 

In this paper, we present a systematic study of strong 
and weak lensing analysis for a sample of 28 clusters from 
the SGAS, based on our extensive fo llow-up imaging obser - 
vations with Subaru/Suprime-cam (|Mivazaki et al.1 12OO2I ) . 
We study the radial dark matter distributions of the clus- 
ters in detail to measure the concentration parameters for 
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these clusters. For this purpose we conduct a stacked lensing 
analysis as well as individual modelling of clusters. With the 
stacked lensing technique we also study the two-dimensional 
mass distribution, and constrain the ellipticity of the pro- 
jected mass distribution. 

The structure of this paper is as follows. We describe 
our cluster sample and follow-up imaging observations with 
Subaru/Suprime-cam in Section [2l and present the strong 
and weak lensing analysis in Section [S] The results are 
combined to discuss the mass-concentration relation in Sec- 
tion m We also conduct a stacked lensing analysis, which 
is detailed in Section [S] We then study the cluster mem- 
ber galaxy distribution in Section [6l and summarize the 
main conclusions in Section [71 In Appendix [A] we conduct 
semi-analytic calculations to evaluate the effect of the lens- 
ing bias on various observables. Appendix [B] summarizes 
the strong and weak lensing analysis results of individual 
clusters. Throughout the paper we assume the standard A- 
dominated flat cosmological model with the matter density 

= 1 — = 0.275, the dimensionless Hubble constant 
h — 0.702, the baryon density Q^hh? — 0.02255, the spectral 
index ris — 0.968, and a normalization for t he matter power 
spectrum as = 0.816 (|Komatsu et al.ll201lh where needed. 



2 CLUSTER SAMPLE AND FOLLOW-UP 
OBSERVATIONS 

2.1 Cluster sample 

We draw our sample of clusters for detailed lensing analysis 
from SGAS, which involved an initial extensive visual search 
of giant arcs in red-sequence selected clusters in the SDSS. 
The sample is constructed utilizing two selection methods. 
One is the SDSS "Visual" survey which selects giant arc can- 
didates from the visual inspection of the SDSS imaging data 
(M. D. Gladders et al., in preparation). The other survey is 
the SDSS "Blind" survey which searches for giant arcs from 
^f-band follow-up imaging of the most massive ^ 200 cluster s 
selected from the SDSS imaging data (|Hennawi et al.ll2008h. 
Some of the clusters have also been reported in IWen etaL 
(l201lh . We are conducting massive spectroscopic follow-up 
observations of these new giant arcs with the Ge mini Multi- 
Object Spectrograph fGMOS: lHook et aLlliooi ) to measure 

redshifts for the newly discovered arcs as well as redshifts 

I 1 II ' 

for the lensing clusters. See Bayliss et al. (2011b) for details 

of our spectroscopic follow-up and successful redshift mea- 
surements for more than 20 giant arcs. We note that this 
paper includes a few new G MOS spectroscopy re sults ob- 
tained after the publication of lBavliss et al.l (|2011bh . We also 
include two cluster-scale quasar lenses discovered from the 
SDSS Quasar Lens Sea rch (SQLS; Oguri et al. 2006, 2008a; 
llnada et al.ll2008l . l2010l ) , as they satisfy similar selection cri- 
teria to the SGAS lens sample, and have also been observed 
at the Subaru telescope. 



2.2 Observations with Subaru/Suprime-cam 

We o bserved the SG AS giant arc clusters with Suprime- 
cam (|Mivazaki et al.l l2002i ) at the Subaru 8.2-meter tele- 
scope, primarily for the wide-field weak lensing analysis, 
between 2007 June and 2011 April. Combined with a few 



imag es retrieved f rom the archive system named SMOKA 
(Bab a et al.ll2002l ). our sample comprises 28 clusters with 
Subaru/Suprime-cam multicolour follow-up imaging. The 
Suprime-cam has a large field-of-view of ^ 34^ x 27' with 
a pixel scale of 0^202, and is therefore ideal for weak lensing 
studies of massive clusters at intermediate redshifts whose 
typical virial radius is < 10^ Our strategy is to obtain deep 
images in r-, and i-bands, with the longest exposure in 
r-band. The deep r-band imaging is performed to conduct 
weak lensing analysis using r-band images, but the addi- 
tional colour information from g- and i-bands are crucial for 
reliable selection of background galaxies for the weak lens- 
ing analysis (more details presented in Section 13. 2p as well 
as secure identifications of multiply imaged arcs in the clus- 
ter cores. We note that the weak lensing analy sis results of 
the first 4 SGAS clusters have been published in I Oguri et al.l 
(2009). 

The data ar e reduced using SD FRED and SDFRED2 
(jYagi et al.ll2002l : lOuchi et al1l2004l l. The reduction proce- 
dure includes bias subtraction, fiat-fielding, distortion cor- 
rection, sky subtraction, and co-adding to generate the fi- 
nal mosaic images. When co-adding, we remove some of the 
frames which have significantly worse seeing sizes than other 
frames. We derive the magnitude zero-points by comparing 
objects in the reduced images with the photometric cata- 
logue of the SDSS data. The astrometric calibration is per- 
formed with SCAMP (^Bertin 2006), again using objects in 
the SDSS data as astrometric reference sources. The pho- 
tometric galaxy catalogue of each cluste r is constructed us- 
ing SExtractor (iBertin k, Arnoutsl 1996lV with the Galacti c 
extinction correction (Schlegel, Finkbeiner, & Davis "l998'). 
Our 28 cluster sample and basic parameters of follow-up 
images are summarized in Table [1] 



3 STRONG AND WEAK LENSING ANALYSIS 
3.1 Strong lensing analysis 

We mostly follow lOguri et al.1 (|2009l l for the strong lens mod- 
elling methodology. Basically we assume the elliptical ex- 
tension of the NEW profile (jNavarro et al.1 Il996l . Il99/t ) to 
model the mass distribution of a dark halo, and add con- 
tributions from member galaxies assuming the pseudo-Jaffe 
model. The velocity dispersion a and cutoff scale rcut of the 
pseudo-Jaffe model are assumed to scale with the luminosity 
as L oc cr^/^ and rcut oc L^^^. The normalization and cut-off 
radius may be fixed to a typical value if there are not enough 
observational constraints. We may also fix the centre of the 
main halo to the position of the brightest cluster galaxy, 
depending on the number of available constraints and the 
configuration of multiple images. 

We identify multiple images based on spectroscopic 
follow-up results of lBavliss et al.l ()2011bh as well as colours 
of lensed arcs measured from the Subaru/Suprime-cam im- 
ages. While redshifts of arcs are available for many of 
our sample clusters, for clusters without any arc spec- 
troscopy information we assume that the redshift of the 
main arc is Za.rc — 2 with the conservative la error of 
Q"(^arc ) — 1, because our extensive spectroscopic observa- 
tions (iBavliss et al]|2011al l as well as photometric analysis 
(|Bavlissll20l'lf )' have convincingly shown that most of the arcs 
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Table 1. Cluster sa mple and summary of Subaru/Suprime-cam imaging observations. The redshifts of the clusters are taken from 
iBavliss et aD (l2011bl V 



Name R.A. Decl. z Exp {g) Seeing {g) Exp (r) Seeing (r) Exp (i) Seeing (i) 

(J2000) (J2000) (sec) (arcsec) (sec) (arcsec) (sec) (arcsec) 
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" Photometric redshifts estimated from the SDSS data, as spectroscopic cluster redshifts are not available for these clusters. 
^ Based on the spectroscopy of the brightest cluster galaxy at Apache Point Observatory 3.5-meter telescope. 

^ We use deep i?c-band images retrieved from SMOKA instead of obtaining r-band follow-up images. This field includes two separate 
cluster cores, both of which act as strong lenses. 

^ Based on the new spectroscopic observation with Gemini/GMOS conducted after the publication of iBavliss et al.l (|2011bl ). 



in our sample fall in this redshift range. The positional un- 
certainties of all the arcs are assumed to be l'' in the image 
plane, which is much larger than measurement uncertainties 
but is typical of the uncertainties associated with parametric 
strong lens modelling of clusters. 

We perform strong lens modelling using the software 
glafic (|Ogurill2010l\ us in^ the m inimization in the source 
plane fseelpgu ri et al.ll2009l : IOguril»2010l From the derived 
best-fit mass model we compute the Einstein radius of 
the system by solving the following equation: 

^(<^e) = 4t/ n{0)dO=l, (1) 

and use only this Einstein radius determination as the con- 
straint from our strong lensing analysis. We use the Einstein 
radius as it is a robust quantity that is well constrained 
from the strong lensing modelli ng and is fairly in sensitive 
to our model assumptions fe.g.. I Julio et alJ 1200 7') . We es- 
timate the error on the Einstein radius by changing one 
(usually mass) of the parameters from the best-fit values, 
optimizing other parameters, and monitoring the differ- 
ences. However, to be conservative in our estimates, we al- 
ways assign a minimum error of 10% to the Einstein radius. 



even if the procedure above returns smaller errors, given the 
unaccounted complexity of cluster mass distributions and 
projections along line-of-sight (e.g., IP'Aloisio NataraianI 
I2OIII ). For clusters without arc spectroscopy, we conserva- 
tively estimate the la error on the Einstein radius by chang- 
ing the arc redshift by =bl. The Einstein radii are computed 
for both the arc redshift ^arc and the fixed source redshift 
of = 2. We use the former value for the combined strong 
and weak lensing analysis, whereas the latter value is used 
for the statistical, stacked lensing analysis. 

We note that the definiti on of the E i nstei n radius 
adopted here differs from that in lOguri et al.1 (|2009l ). In this 
new definition th e contribution of st ellar masses is explicitly 
included, while in lOguri et al.1 (|2009l ) the Einstein radius was 
computed from the dark halo component alone. The use of 
this modified definition is because the separation of dark 
and luminous components is not obvious, and has to rely on 
several assumptions, unless we have additional data such as 
the velocity dispersion measurement of the central galaxy 
(e.g.. San d et al.1 l2008l : INewman et al.ll2009l ). Thus the ef- 
fect of baryonic components has been taken into account in 
interpreting our results. 

Table [2] lists the Einstein radii derived from strong 
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Table 2. Summary of strong lensing analysis. We tabulate the Einstein radii both for arc redshifts and for = 2, the total number of 
multiple images used for strong lens modelling (A^img; the number of multiple image sets is shown in parentheses), the best-fit ellipticity 
(e) and position angle (^e) of the main dark halo component, and references for arc redshifts. 



Name z^rc"" A^img 6'e {zs = ^arc) 6'e {zs =2) e Oe^ Ref.^ 

(arcsec) (arcsec) (deg) 



SDSSJ0851+3331 


1.693 


4(1) 


21.6+?-? 

— 2.2 


23.0+^-^ 

— 2.3 





23 


39.1 


1 


SDSSJ0915+3826 


1.501 


3(1) 


— 1.0 


11.4+;? 

— 1.1 





28 


-74.9 


1,2 


SDSSJ0957+0509 


1.820 


3(1) 


5.2+0-^ 

— 0.5 


5.4+°1 

— 0.5 





82 


64.2 


1,3 


SDSSJ1004+4112 


1.734 


31(8) 


— 0.7 


— 0.9 





14 


-28.0 


4,5 


SDSSJ1029+2623 


2.197 


3(1) 




9.9+tl 





58 


-81.5 


6,7 


SDSSJ1038+4849 


2.198 


9(3) 


12.6+}r 

— 1.6 


n.2+ii 

— 1.4 





15 


-52.8 


1,8,9 


SDSSJ1050+0017 


2±1 


3(1) 




16.1+i^,° 

— 2.5 





28 


-6.9 




RCS2J1055+5548 


1.250 


3(1) 


10.0+}° 

— 1.0 


12.7+1-^ 

— 1.3 





82 


3.2 


1 


SDSSJlllO+6459 


2±1 


3(1) 


8.4+ 

— 1.5 


8.4+]^'' 

— 1.5 





60 


76.9 




SDSSJ1115+5319 


2±1 


5(1) 


21.9+^^.^ 

— 3.4 


21.9+^^.^ 

— 3.4 





63 


-47.4 




SDSSJ1138+2754 


1.334 


5(2) 


9.8+;° 

— 1.5 


12. 9+^-^ 

— 3.6 





50 


-18.7 


1 


SDSSJ1152+3313 


2.491 


7(2) 


8.7+°-9 

— 0.9 


o 9 + 0.8 
— 0.8 





34 


64.5 


1 


SDSSJ1152+0930 


2±1 


3(1) 


— 0.7 


4.5+8-2 

— 0.7 





80 


-38.5 




SDSSJ1209+2640 


1.021 


6(2) 


8 8+°-^ 

— 0.9 


2I.3 + 0I 
— 2. 1 





14 


66.1 


1,10 


SDSSJ1226+2149 


1.605 


3(1) 


— 2.6 


15.7+?-« 

— 2.8 





47 


-57.5 


1 


SDSSJ1226+2152 


2.923 


3(1) 


10.0^2-^ 







18 


12.3 


1,11 


A1703 


2.627 


21(6) 




25.3+25 





33 


-26.2 


1,12,13 


SDSSJ1315+5439 


2±1 


3(1) 


16.9+^9^ 


16.9+^9^ 





18 


-36.8 




GHO132029+3155 


2±1 


4(1) 


21.5+^-0 


21.5+^-0 





10 


76.5 




SDSSJ1329+2243 


2.040 


3(1) 


io-9t}:6 


io-8t}:6 





22 


-4.9 


14 


SDSSJ1343+4155 


2.091 


3(1) 


5 4+2.5 







73 


64.9 


1,2,15 


SDSSJ1420+3955 


2.161 


6(2) 


Q Q + 2-8 
^•^-1.0 


9 4+2.7 





77 


73.6 


1 


SDSSJ1446+3032 


2±1 


4(1) 




16.8 + ^33! 





22 


-41.6 




SDSSJ1456+5702 


0.833 


6(2) 


13.2+}-^ 


30 1+^-2 





21 


-13.6 


1 


SDSSJ1531+3414 


1.096 


6(2) 


11 7+^-2 


16.6^}-; 





09 


-37.5 


1 


SDSSJ1621+0607 


4.131 


3(1) 


12-5+25 


10-3^22 





33 


-40.0 


1 


SDSSJ1632+3500 


2±1 


4(1) 




14.3+^Y 





30 


-61.9 




SDSSJ2111-0114 


2.858 


9(3) 


-17 7+10.6 
'-5.8 


14 5+^-2 





24 


10.0 


1 



" Note that for some clusters there are multiple arcs with different source redshifts. In this case we show the redshift of the main arc. 
^ The position angle Oe is measured East of North. 

^ 1 -iBavliss et al.l (1201 Ibl): 2 -iBayliss et al.l (l2010h: 3 -iKubo et al.l (l2010h : 4 -llnada et al.l (l2003h : 5 -ISharon et al.l (|2005|): 6 - 
Inada et al.l (|2006|); 7~ ]Og^uri et al.l (l2Q08bl): 8 iBelokurov et al.l (l2009l); 9 -iKubo et al.l (120091 )710 - lOfek. Seitz. KleinI J2008r) : 11 - 

new redshift not yet published in the literature; 15 



inada et al.l (l2UUbl): 7 - lUguri et al.l (\200m): b -IJjelokurov et al.l (l2U(jyi): 9 -IKub 
Koester et ai.l (I2010I ): 12 iLimousin et al.l (20081) : 13 - iRichard et al.l (|2009l ): 14 



Diehl et al.l (|200' 



lensing modelling. SDSSJ1226+2149 and SDSSJ1226+2152 
are two nearby strong lensing cores in a highly com- 
plex massive structure, both of which are covered by our 
Subaru/Suprime-cam imaging (Bavlis s et alj 1201 Ibl V We 
perform strong lensing analysis for both cores. We also show 
the best-fit ellipticity e and the position angle Oe of the 
best-fit main NFW component, as they are used later for 
two-dimensional stacking analysis fSection l5.2|) . The best-fit 
critical curves for individual clusters are presented in Ap- 
pendix [B] 

Strong lensing analysis for some of th ese clusters 
has been published in t he literature (e.^. , lOguri et alJ 
'2009; Koester et alj |201Q| : 'Bavliss et all I2OIOI : lOguril |201Q| : 
,Gralla et al. 201 J). Our estimates of the Einstein radii pre- 
sented here generally agree with these previous results. For 
A 1703, our best-fit Einstein radius is in excellent agree- 
me nt with that of the independent strong lens modelling 
bv lRichard eTaD (|20Q9l ). An exception is SDSSJ1343+4155 



whose Einstein radius in our strong le ns modelling r e sult is 
much sm aller tha n what presented in 'Bavliss et al.l (|201Qh 
and Gra lla et al.l (|2011) because a large offset of the halo 
centre from the brightest cluster galaxy has been found in 
the previous analysis, whereas in the present paper the cen- 
tre of the halo is fixed to the location of the brightest cluster 
galaxy. 

3.2 Weak lensing analysis 

For weak le nsing measur ements, we follow the formalism 
outlined in iKaiser et all (KSB; Il995l ) using the software 
package IMCAtQ. We first detect objects in the reduced 
images using a hierarchical peak finding algorithm. For all 
the clusters in our sample, we use r-band images for weak 

^ http://www.ifa.hawaii.edu/ kaiser/imcat/ 
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lensing shear measurements. We then measure the shapes 
of objects by iteratively refining the centroid of each ob- 
ject. Stars for correcting the distortion of the Point Spread 
Function (PSF) are selected in a standard way by identify- 
ing the appropriate branch in the magnitude-half light ra- 
dius Th plane, along with the peak significance cut > 15. 
Shapes of the stars are measured as a function of the size 
of the weight function, r^, in order to make PSF correc- 
tions with matched Vg values. We divide each co-added im- 
age into 4x3 chunks and fit the PSF in each chunk in- 
dependently with second order bi-polynomials. The smear 
polarizability is corrected by computing a scalar polariz- 
ability Ps from the trace of the matrix, and then fitting 
Ps as a function of magnitude, r^, and the galaxy ellip- 
ticity. For our weak lensing analysis, we only use galaxies 
with f > lb and r^, > + 2cr(r^), where and cr(rf^) are 
median and root-mean-square dispersion of half-light radii 
for the stars selected above. Given the general tendency 
of the KSB algorithm t o underestimate the weak lensing 



shear (|Erben et al.ll200ll : lHevmans et al.ll2Q06l : lMassev et al 
l2007l ). we include a calibration factor of 1/0.9, i.e., this fac- 
tor is multiplied to all the estimated shear values. For each 
object we assign the statistical weight Wg defined by (e.g., 
Hamana et al. "2003"; 'Mivaza ki et al.l l2007l: [Hamana et al. l 
200 9: Okabe et al. 2010; Oka be. Okura. Futamasd l2010l : 
Umetsu et al.ll201Q . l2011b ). 



1 



(2) 



with a = 0.4 and ag is the variance of the shear computed 
from 20 neighbors in the magnitude-r^ plane. When comput- 
ing shear by averaging shear measurements of galaxies in a 
bin, we also apply a 3a clipping which appears to reduce the 
shear measurement bias. 

To check the accuracy of the weak lensing shear mea- 
surement, we perform a series of image simulations. Specif- 
ically, we generate a galaxy catalogue using the software 
Stuff (|Bertinll2009l ). Each galaxy is described by the sum 
of bulge and disk components, which we model with Ser- 
sic profiles with the index n = 4 and n = 1, respectively. 
We also add stars in the catalogue. We convolve the im- 
age with a PSF which we assume follows the Moffat pro- 
file S(r) oc [1 + (r/a)^]~^ with an elliptical extension. 
Based on the catalogue, we generate a number of realis- 
tic Subaru/Suprime-cam like images with different seeing 
sizes ( QV5-1'.'1) an d /3 (3 < /3 < 12) using the software 
glafic (jOguril Eoid V We find that th e resulting shear mul- 
tiplicative error (th e parameter m in iHevmans eralll2006l : 
iMassev et al .' 2007) depends on both seeing size and /3 such 
that m is smaller for larger seeing sizes or smaller (3, but 
our algorithm generally yields \m\ < 0.05 for a wide range 
of PSF parameters examined here. 



3.3 Galaxy selection in colour-colour space 

A careful selection of background galaxies is essential for 
cluster weak lensing studies, because contamination by clus- 
ter member galaxies is known to dilute the detected weak 
lensing sig nal significantly, particu l arly near the cluster cen - 
tres (e.g., iBroadhurst et al'l l2005l : iMedezinski et al.ll200/t ). 
Our ^ri-band imaging is very powerful for reliable back- 



I 
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all 
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Figure 1. Upper: The colour cut selecting background galaxies 
for weak lensing analysis (shading). The colour cut for red mem- 
ber galaxy selection is indicated by dotted lines. The four regions 
correspond to colour cuts for different cluster redshifts. Solid con- 
tours indicate galaxy number density in the COSMOS catalogue 
(i < 25). lower: Photometric redshift distributions of galaxies 
in the COSMOS catalogue, before (dotted) and after (solid) the 
colour cut. 



ground galaxy selection, because we can select galaxies effi- 
ciently in colour-colour space (Medezinski et al. 2010). 

We determine the colour cut that is appropriate for our 
cluster sa mple based on the COSMOS photometric galaxy 
catalogue (jllbert et al.l l2009l V Thanks to the wide wave- 
length coverage from the ultraviolet to mid-infrared, the 
photometric redshifts are very accurate down to i ^ 25. 
By inspecting the photometric redshift distributions in each 
point of the g — r versus r — i colour space, we determine 
the colour cut for our analysis as 



g-r>l r-i< OA(g - r) - 0.5, 
g-r < 0.3, 
r-i> 1.3, 
r — i > g — 



(3) 
(4) 
(5) 
(6) 
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Figure 2. The radial number density distributions of galaxies, 
which are obtained by averaging distributions over all 28 clus- 
ters. We show distributions for background galaxies (open cir- 
cles), cluster member galaxies (filled squares), and "field" galaxies 
(filled triangles) that are simply defined as being neither member 
nor background galaxies. Note that the distributions are normal- 
ized by the overall number densities. 



Figure [T] shows the cut and resulting COSMOS photometric 
redshift distribution. It is seen that our cut efficiently selects 
galaxies at ^ > 0.7, the redshifts higher than any cluster 
redshift s in our sample. In our weak lensing analysis, we also 
limit the range of r-band magnitude to 21 < r < rum, where 
the limiting magnitude rum is determined from the galaxy 
number counts of each cluster field image (see below). 

In addition to background galaxies, we identify cluster 
(red) member galaxies by the following criteria: 



ai 



0.5( 



0.3 < ^-r < ai + 0.3, 



(7) 



+ as - 0.15 <r-i< 0.5(g - r) + as + 0.15, (8) 



with (ai, as) are (1.5, -0.25) for z < 0.35, (1.6, -0.22) for 
0.35 <z < 0.45, (1.6, 0) for 0.45 <z< 0.55, and (1.7, 0.05) 
for 0.55 < z. These cuts are also indicated by the upper 
panel of Figure (TJ 

Figure [2] shows the average radial number density dis- 
tributions of galaxies, background, member, and field (i.e., 
galaxies which are neither background nor member) galaxies 
in our 28 cluster sample. As expected, the member galaxy 
density is steeply rising near the cluster. Field galaxies do 
not show strong dependence on the distance from the centre, 
as that should be dominated by local foreground galaxies. A 
slight increase at small radius suggests that some blue clus- 
ter member galaxies are included in the field galaxy sample. 
On the other hand, background galaxy density decreases 
near the centre, which is in fact expected because of the 
lensing magnification and dilution effect of clusters (e.g., 
iBroadhurst et al.ll2005l : lUmetsu et~alll2011al ). These distri- 
butions basically support the successful colour cut for se- 
lecting background galaxies. 

In order to extract physical quantities for clusters from 
weak lensing signals, we need an estimate of the lensing 
depth of our source galaxy sample. Again, the COSMOS 
photometric redshift catalogue is useful for this purpose, but 
one problem is that many of our source galaxies are fainter 




Figure 3. The mean galaxy redshift as a function of i-band mag- 
nitude, which is derived from the COSMOS photometric redshift 
catalogue ( Ilber t et al...2009i V The solid line plots our fit given by 
equation ()10|) . 



than the magnitude limit of the COSMOS photometric red- 
shift catalogue. To overcome this problem, we adopt the fol- 
lowing procedure to estimate the lensing depth . First, from 
the COSMOS photometric redshift catalogue (jllbert et al.l 
2009) we derive the response function for our colour cut 
as (l)cut(z) = ncnt(z)/ntot(z), where ntot(z) and ncut(z) are 
photometric redshift distributions of galaxies before and af- 
ter the colour cut, respectively. Our weak lensing analysis 
involves the statistical weight Wg (equation [2]) , and hence 
we compute the mean weight as a function of i-band magni- 
tude as well, w^(i), from the background galaxy catalogue. 
Next we derive the redshift distribution of galaxies as a func- 
tion of the i-band magnitude as follows. Fir st, we adopt the 
functional form proposed bv lSchrabback et al.1 (|2010|): 



P{Z\Z) oc 



+ cu exp 



(- 



,(9) 



with u = max(0,'i — 23) and (a, (3, c, d, 7) = (0.678, 5.606, 
0.581, 1.851, 1.464). We then recompute the mean galaxy 
redshift as a function of i-band magnitude for the range 
17 < i < 25 using the COSMOS photometric redshift cata- 
logue. The result shown in Figure [3] suggests that it is well 
approximated by 



0.1(i- 
0.25(i 



22) +0.6 
-22) + 0.6 



(17<i<22), 
(22 < i), 



(10) 



We extrapolate this linear relation to a fai nter magnitude 
of i > 25, because ISchrabback et al.] (|2010h has explicitly 
shown that such extrapolation can explain the redshift dis- 
tribution of galaxies reasonably well down to i ^ 27. Us- 
ing this relation, we derive the relation between zo in equa- 
tion Q and i-band magnitude, which is approximated by 



( 0.16(i - 22) + 0.97 (17 <i<22), 

OA(i - 22) + 0.97 (22 <i< 23.1), 

1.41 (23.1 < i < 24.3), 

^ 0.2(i - 24.3) + 1.41 (24.3 < i). 



(11) 



Using this redshift distribution, we derive the mean depth 
as 
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' Dis 



,.dN 



(12) 



Di 



I 



dz^p{z\i)(l)c 

J-^os 



J dzp 



{Z\i)^cut{z) ] ,(13) 



where dN/di are the z-band number counts of background 
galaxies used for weak lensing analysis. We define the effec- 
tive source redshift ^s,eff such that it reproduces the mean 
depth: 



Dis . . 

J-^os 



/ Dis 

\Dos 



(14) 



Throughout this paper, we assume that all the galaxies are 
located at ^s,eff for our weak lensing analysis. 

The depth estimate above can be affected by the lensing 
magnification, because the magnification enhances the effec- 
tive lensing depth. However, the magnification factor tends 
to decrease rapidly beyond the Einstein radius such that a 
typical magnification factor of the innermost radial bin for 
our weak lensing analysis is ^ 20% or so. This corresponds 
to the enhancement of the shear amplitude of ^ 5%, which 
is not significant compared with other uncertainties. More- 
over the effect is much smaller at radii where weak lensing 
signals mostly come from. Thus we conclude that the mag- 
nification effect is insignificant for our analysis. Photometric 
errors that become larger for fainter galaxies can smear the 
distribution in the color-color space, and hence can affect the 
lensing depth estimate, but photometric errors are included 
in the COSMOS photometric redshift catalog as well. Pho- 
tometric errors for the faintest galaxies in our weak lensing 
analysis are similar to those for the faintest galaxies in the 
COSMOS catalog. Another possible systematic effect is the 
size cut, which may systematically eliminate galaxies at high 
redshifts and hence bias the lensing depth estimate. While 
we expect the effect to be small because of the conservative 
magnitude limit we adopt and the narrow redshift distribu- 
tion after the color cut, we leave the detailed exploration of 
this effect for future work. 

Table [3] summarizes our weak lensing analysis. It is 
seen that our PSF correction algorithm successfully reduces 
the stellar ellipticities to the level that are much smaller 
than typical lensing shear amplitudes for our cluster sam- 
ple (7 ^ 0.01). The relatively small number density of 
Ug ~ 10 arcmin"^ is due to the colour cut for selecting 
background galaxies. From the shear catalogues we com- 
pute tangential shear profiles. The tangential shear is 
computed from the reduced shear g = (gi, ^2) as 

^+ = — ^1 cos 20 — ^2 sin 20, (15) 

where is the polar angle. Throughout the paper, we assume 
that the position of the brightest galaxy in the strong lensing 
region is coincident with the cluster centre. Although our 
strong lensing analysis suggests that the position of the mass 
peak can differ slightly from the central galaxy position, we 
find that any such offset is much smaller than the typical 
inner radial boundary of our weak lensing analysis. Hence 
the effect of off-centreing on our results should be negligibly 
small. The average shear value in each radial bin is estimated 
by the weighted mean of the tangential shear as follows: 



9+ 




(16) 



where the summation runs over all galaxies in the bin. Sim- 
ilarly the statistical error in the tangential shear measure- 
ment in each bin is computed from tlie__weiglited average of 
the variance of the shear, cr^ (e.g.JOkabeetal. 2010). The 
total signal-to-noise ratio of the tangential shear profile is 
shown in Table [3] and indicates that weak lensing signals 
are detected significantly (S/N > 5) for most of our clus- 
ters. The mass maps reconstructed from the weak lensing 
analysis are shown in Appendix [B] 



4 COMBINING STRONG AND WEAK 
LENSING 

4.1 Methodology 

We combine constraints from the tangential shear profile 
with constraints from strong lensing. As in I Oguri et al.l 



(120091 ). we combine both sets of constraints by summing 
up x^: 

2 2,2 

X =XSL + XWL- 



(17) 



We include strong lensing constraints from the Einstein ra- 
dius at the arc redshift: 



2 _ [fe -6'E(Mvir, Cvir)]^ 
XSL - -o 



(18) 



where and cte are the best-fit Einstein radius at Zs = ^arc 
and its error presented in Tabled and fe(Mvir, Cvir) is the 
predicted Einstein radius assuming the NEW profile. On 
the other hand, weak lensing constraints come from binned 
tangential shear measurements: 



2 [9+,i - 9+{0i'i Mvir, Cvir)]^ 



XWL 



(19) 



where g^^i and ai are observed reduced shear and its error 
at i-th radial bin, and —g^iOi] Mvir, Cvir) is the predicted 
reduced shear by the NEW model. The NEW profile, which 
we adopt as an analytical model for the radial mass distri- 
bution, has the three-dimensional density profile of 

Ps 



p{r) = 



where 



ps 



(r/r.)(l + r/r.)2' 



Avir(2;)pm(^)c^ir 



(20) 



(21) 



3 [ln(l + Cvir) - Cvir/(1 + Cvir)] ' 

where A(z) is the nonlinear over-density predicted by the 
spherical collapse model. We parametrize the profile with 
two parameters, the virial mass Mvir 



Mvir = ^rl;^/\viv{z)pm{z), 

and the concentration parameter Cv 

Tvir 



(22) 



(23) 



4.2 Correlation between strong and weak lensing 
observables 

Before deriving best-fit parameters for individual clusters, 
we first compare strong and weak lensing observables. 
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Table 3. Summary of the weak lensing analysis. We show the median of stellar ellipticities before (e*^^) and after (e*oj.) the PSF 
correction (numbers in parentheses are standard deviations), the number of stars used for the PSF correction (N*), the number density 
of source background galaxies (rig), the r-band magnitude limit of source background galaxies (nim), the mean lensing depth ({Dig/ Dos)) 
and the effective source redshift {zs^efi)^ the range of the radii for tangential shear fitting (0 range), and the total signal-to- noise ratio of 
the tangential shear profile (S/N). 



Name e*^^ e*^^ A^* rig (Dis/Dos) Zs,ef¥ range S/N 







n n 




) 






— 


) 




(arcmin ^) 


(m3 












(arcmin) 




SDSSJ0851+3331 


1 


51 





90) 





01 




47) 


834 


10.1 


25 


8 





595 


1 


098 


0.63- 


-15.85 


6.93 


SDSSJ0915+3826 


1 


40 





83) 





02 


;o 


54) 


694 


11.7 


25 


6 





571 


1 


109 


0.40 


-6.31 


4.59 


SDSSJ0957+0509 


1 


28 


1 


74) 





01 


;o 


93) 
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namely the Einstein radius and the virial mass, which should 
serve as a sanity check for results of more detailed analysis 
presented in the following sections. For the strong lensing 
observable, we use the Einstein radius for the fixed source 
redshift of Zs = 2 in order to compare results for different 
clusters. Given the strong degeneracy between Mvir and Cvir 
in weak lensing analysis, we quantify the strength of the 
weak lensing signal by fitting the tangential shear profile 
for a fixed concentration parameter of Cvir = 6, a typical 
value for the concentration when the lensing selection effect 
is taken into account (see below) . Figure [H shows the corre- 
lation between the Einstein radius from strong lensing and 
the virial mass from weak lensing. We confirm that these two 
measurements are indeed correlated with each other such 
that the more massive clusters detected via their weak lens- 
ing have larger Einstein radii on average. 

Assuming the NEW profile, the relation between mass 
and Einstein radius is determined by the concentration pa- 
rameter. We also plot the expected correlation for the con- 
centration assuming Cvir = 6, as well as the scatter of the re- 
lation originating from the shift of the concentration by =b2, 
for the lens redshift z = 0.46 which corresponds to the me- 
dian redshift of our cluster sample. We find that the observed 
correlation roughly follows the expectation for Cvir = 6, al- 
though there is a tendency for the Einstein radii to be larger 
than expected, particularly for low-mass clusters. This anal- 



ysis suggests that our sample of clusters are slightly more 
concentrated than Cvir 6, and that clusters with lower 
masses on average have higher concentration values. 

4.3 Fitting results 

Tabled] summarizes the results of the combined analysis, i.e., 
simultaneous fitting of the observed Einstein radius and tan- 
gential shear profile with the NEW model predictions. The 
tangential shear profiles for individual clusters are compared 
with best-fit models in Appendix [B] Fitting is performed in 
the parameter range IO^^H'^Mq < Mvir < lO^^h'^M^ and 
0.01 < Cvir < 39.8 1. For four clusters studied previously in 
lOguri et al.1 (|2009l \ our new results are fully consistent with 
the old result. In addition, the best-fit mass and concentra- 
tion parameter for A 1703 are in good agreement with the 
result of an independent strong and weak lensing analysis 
bv lZitrin et al.1 (|201Q| ). 

With a large sample of clusters with measurements 
of the concentration parameter from combined strong and 
weak lensing, we can study the mass-concentration relation 
quite well. A caveat is that our sample of clusters are se- 
lected as those having prominent arcs. It has been noted 
that clusters selected by strong lensing (e.g., by giant arcs 
or large Einstein radii) represent a strongly biased popu- 
lation such that the concentration parameter inferred from 
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Table 4. Constraints on the mass Mvir and concentration parameter Cvir from weak lensing and combined strong and weak lensing 
analysis. Errors indicate Icr errors on each parameter after marginahzing over the other parameter. 
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Figure 4. Correlation between strong and weak lensing observ- 
ables. The best-fit mass Mvir for a fixed concentration of Cvir = 6 
from just weak lensing is compared with the Einstein radius for 
the source redshift Zs = 2 from strong lensing. The red line with 
shading shows the expected M^i-^-Oe relation assuming the NFW 
profile with the concentration of Cvir = 6 zb 2. 



the projected mass distribution is on average much larger 
mostly due to the halo triaxialitv (IHennawi et ah! l2007l : 
iQguri &; Blandfordll2009l : iMeneghetti et al.ll2010l ) .Although 
m 'Qg uri et alJ (|2009l " we assumed the constant enhancement 
of the concentration parameter due to the lensing bias, sim- 
ple considerations suggest that the lensing bias of the mass- 
concentration relation should depend strongly on the mass. 
To derive more accurate theoretical predictions based on the 
ACDM model, in Appendix [A] we conduct a series of semi- 
analytic calculations with ray-tracing of extended sources 
to estimate the effect of the l e nsing bias, based on a triax- 
ial halo model of IJing Sutol (|2002i). For reference, we find 
that the mean mass-concentration relation at z = 0.45 with 
the lensing bias predicted by this model is roughly described 
by 



cvir(^ = 0.45) ^ 6.3 



Mv 



5 X 10i4/i-iM( 



(24) 



which show relatively strong dependence on the halo mass, 
simply because of the mass dependence of the lensing bias 
(see Appendix [A}. The scatter of the relation is estimated 



to be (Tio 



: 0.12. 



Figure [5] shows the mass- concentration relation ob- 
tained from our lensing analysis for 28 systems. For com- 
parison, we also show accurate lensing measurements of the 



© RAS, MNRAS 000,[TH22] 



Combined strong and weak lensing analysis of 28 clusters 11 



20 - 



— 'k II f^TTl 1 r 



-1 I 




12007 



.Buote et aD l2007l : lEttori et all \201(t ) analysis (see 
Bivianollioosl ). Our result suggests that the observed 



1 10 
M,,, [101%-iMj 

Figure 5. The mass-concentration relation obtained from com- 
bined strong and weak lensing analysis. Filled triangles show our 
results presented in this paper, whereas filled squar es show results 
from l iterature; A1689 , A370, CL0024 RXJ1347 (jUmetsu et al.l 
l2Qllbr) . and A383 (Zi trin et all l2Qllb h. The black shaded re- 
gion indicates the predicted concentration parameters as a func- 
tion of the halo mass with the lensing bias taken into account 
(see Appendix |A] for details). The solid line is the best-fit mass- 
concentration relation from fitting of our cluster sample (i.e., filled 
triangles), with the la range indicated by dotted lines. 



concentrati on parameters for 5 massive clusters fro m the 
literature (jUmetsu et alJ l2011bl : IZitrin et al.1 l2011bl ). The 
Figure clearly indicates that measured concentrations are 
correlated well with the mass. More massive clusters have 
on average smaller concentrations, which is consistent with 
the theoretical expectation, although the slope is obviously 
much steeper. 

To put this on a more quantitative footing, we fit 
the mass-concentration parameter with both the normal- 
ization and the mass slope as free parameters. Here we ig- 
nore the redshift dependence of the mass-concentration re- 
lation, given the predicted little evolution of the concentra- 
tion of massive ha loes with redshift fe.g. JZhao et al.ll2009l : 
IPrada et al.ll201l1 V We use the following estimator for fit- 
ting: 



X 



: (Cvir,obs) - log (Cvir,fit)]^ 



(25) 



where cTgt is the la measurement error on log(cvir) for indi- 
vidual clusters from the lensing analysis (see Table (D) , and 
(Jin = 0.12 is the intrinsic dispersion predicted by our calcu- 
lations. We find that the best-fit mass-concentration relation 
from our lensing sample of 28 clusters is 

cvir = (7.7 ± 0.6) , (26) 



^5 X lOi^/i-iMoy/ 

where we also included la errors on the normalization and 
the slope. 

We detect a strong mass dependence of the concentra- 
tion parameter with the slope of —0.59 =b 0.12 in our lensing 
sample, which should be compared with the predicted slope 
of ^ —0.2 for the strong lensing selected sample of clusters. 
We note that steeper mass-concentration relations than the- 
oretical expectations were a lso suggested by previous weak 
lensing (|Okabe et al.1 l2010l ) and X-ray (jSchmidt &: AllenI 



also 

mass-concentration relation is in reasonable agreement with 
the simulation results for very massive haloes of Mvir 
10^^h~^M(7). The agreement may be ev en better if w e adopt 
recent results of A/'-body simulations bv lPrada et al. I (|201lh . 
who argued that previous simulation work underestimated 
the mean concentrations at high mass end (see also Ap- 
pendix|A}. In contrast, we find that observed concentrations 
are much higher than theoretical expectations for less mas- 
sive haloes of Mvir - IO^^H-^Mq, even if we take account 
of the mass dependence of the lensing bias. 

A possible concern is the correlation of Mvir and Cvir for 
fitting of individual clusters which has been ignored in deriv- 
ing the mass-concentration relation. We examine the possi- 
ble effect of the degeneracy between Mvir and Cvir by the fol- 
lowing Monte Carlo simulation. For each cluster, we change 
best-fit values of Mvir and Cvir by randomly picking up a 
point within the la confidence region in the Mvir-Cvir plane. 
After choosing new best-fit parameters for all the clusters, 
we re-fit the mass-concentration relation to derive best-fit 
values of the normalization and slope in equation (|26p . We 
repeat this simulation for 300 times to check how the corre- 
lated errors between Mvir and Cvir can affect our conclusion. 
We find that the resulting distribution of the best-fit normal- 
ization value is 7.6=b0.2, and that of the slope is — 0.56=b0.05, 
which are small compared with the statistical errors shown 
in equation (|26|) , suggesting that the effect of the corre- 
lated error is not very significant. We note that lOkabe et al.1 
(|2010h also explored potential impacts of the degeneracy 
between Mvir and Cvir on their mass-concentration measure- 
ment very carefully, and concluded that the effect is insignif- 
icant, which is consistent with our finding. 

There are a few possible explanations for the excess 
concentration for small mass clusters. Perhaps the most sig- 
nificant effect is baryon cooling. The formation of the central 
galaxy, and the accompanying adiabatic contraction of dark 
matter distribution, enhances the core density of the clus- 
ter and increases the concentration parameter value for the 
total mass distribution. This effect is expected to be mass 
dependent such that lower mass haloes are affected more 
pronouncedly, simply because the fraction of the mass of 
the central galaxy to the total mass is larger for smaller 
halo masses. Indeed, simulations with radiative cooling and 
star formation indicate that the concentration can be signifi- 
cantly enhan ced by baryon physics particularly for low- mass 
haloe s (e.g., iRudd, Zentner, Kravtsovl l2008l : iMead et al.l 
l2010h , altho ugh the eflFect strongly depends on the effi ciency 
of feedback (Duffy et al. I I2OIOI : iMcCarthv et al. II2OIOV Thus 
baryon cooling appears to be able to explain the observed 
strong mass dependence at least qualitatively. More quanti- 
tative estimates of this effect need to be made using a large 
sample of simulated clusters with the baryon physics as well 
as the proper feedback model included. 



5 STACKING ANALYSIS 

5.1 Stacked tangential shear profile 

We can study the average properties of a given sample by 
stacking lensing signals. This stacked lensing analysis has 
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Table 5. Summary of stacked tangential shear analysis 



0.1 1 
r [h~^Mpc] 

Figure 6. The stacked tangential shear profile obtained by com- 
bining the 25 clusters. The average differential surface density 
{AS+(r)} (see equation I27|) is plotted as a function of the phys- 
ical radius r. Grey points indicate stacked tangential shear mea- 
surements from weak lensing that are not used for fitting. The 
upper left point with a horizontal error-bar is the constraint from 
the average Einstein radius. The solid line with shading is the 
best-fit NFW model with Icr error range. The lower panel plots 
the stacked profile of the 45° rotated component, (AEx(r*)}. 



been successful for constraining mean dark matter distri- 

2006bl: 



butions of cluster samples (e.g;, 
Johnston et~all l2007l : |l 



Mandelbaum et al 

Eli 



Okabe et alJ 



.Leauthaud et alJ 

20101 1. Here we conduct stacking analysis of the tan- 
gential shear profile for our lensing sample for study- 
ing the mass-concentration relation from another view- 
point. Note that the off-centreing effect, which has been 
known to be one of the most signific ant systematic 
errors in stacked lensing an alysis (e.g., Johnston et al. 



2007 : lMandelbaum. Seliak. H irata 2008; Oguri & Takada 
201 ih , should be negligible for our analysis, because of the 
detection of weak lensing signals for individual clusters and 
the presence of giant arcs which assure that selected cen- 
tres (positions of the brightest galaxies in the strong lensing 
region) indeed correspond to that of the mass distribution. 

We perform stacking in the physical length scale. Specif- 
ically, we compute the differential surface density AS+(r) 
which is define by 



AS+(r) = Scr^+((9 = r/L>or 



(27) 



where Scr is the critical surface mass density for lens- 
ing. We stack AS+(r) for different clusters to obtain the 
average differential surface density. We do not include 
SDSSJ1226+2149 and SDSSJ1226+2152 in our stacking 
analysis, because these fields clearly have complicated mass 
distributions with two strong lensing cores separated by only 
- 3'. Furthermore, we exclude SDSSJl 110+6459 as well be- 
cause the two-dimensional weak lensing map suggests the 
presence of a very complicated mass distribution around the 
system. We use the remaining 25 clusters for our stacked 
lensing analysis. 

It should be noted that the reduced shear ^+ has a non- 
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Unear dependence on the mass profile. In fact, the reduced 
shear is defined by ^+ = 7+/(l — /^), where 7+ and n are 
tangential shear and convergence. Thus, the quantity defined 
by equation (|27)) still depends slightly on the source redshift 
via the factor 1/(1 — At), particularly near the halo centre. 
Thus, in comparison with the NFW predictions, we assume 
the source redshift of Zs — 1.1, which is the typical effective 
source redshift for our weak lensing analysis (see Table [3]) . 
Also the non-linear dependence makes it somewhat difficult 
to interpret the average profile, and hence our stacked tan- 
gential profile measurement near the centre should be taken 
with caution. 

It is known that the matter fluctuations along the line- 



of-sight contributes to the to t al e rror budget (e .g., tH oekstra 
2003 : iHoekstra" et al.1 I2OIII : iDodelson 200i; ICmen et al 



201 it ). While we have ignored this effect for the analysis of 
individual clusters presented in Section 2] here we take into 
account the error from the large scale structure in fltting 
the stacked tangential shear profile by including th e full co- 
variance between differen t radial bins (see Oguri Takadal 
I2OIII : lUmetsu et al.ll2011bl , for the calculation of the covari- 
ance matrix). We, however, comment that the error of the 
large scale structure is subdominant in our analysis, because 
of the relatively small number density of background galax- 
ies after the colour cut (see also Oguri et al. 2010). 

In addition to weak lensing, we stack strong lensing 
constraints simply by averaging the Einstein radii for the 
fixed source redshift Zs — 2. This constraint is combined 
with the stacked tangential shear profile from weak lens- 
ing to obtain constraints on the mass and concentration 
parameter for the stacked profile. Note that the Einstein 
radius is related with the reduced shear as ^+(^e) = 1. 
Given the uncertainty from the non-linearity of the reduced 
shear and the the possible bias coming from the uncer- 
tainty of the outer mas s profile (jOguri Hamanal [20T1I : 
iBecker Kravtsovl [20TTI ) . we restrict tangential shear fit- 
ting in the range 0.158/i~^Mpc < r < 3.16/i~^Mpc. However 
we note that our results are not largely changed even if we 
conduct fitting in the whole radius range. 

Figure [6] shows the stacking result for all the 25 clusters. 
The mean cluster redshift for this sample is {z) — 0.469. 
The total signal-to-noise ratio in the whole radius range 
of 0.063/i"^Mpc < r < 5.01/i"^Mpc is S/N = 32. We 
find that stacked tangential shear profile from weak lens- 
ing is fitted well by the NFW profile over a wide range 
in radius. The average Einstein radius from strong lens- 
ing ((6'e) = is slightly larger than the best- 
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Figure 7. Similar to Figure [21 but the stacked lensing analysis 
in three bins is presented. From top to bottom, results for 
largest to smallest bins are shown. Curves and points for the 
largest and smallest ^e bins are shifted vertically by ±0.5 dex 
respectively for illustrative purposes. 
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Figure 8. Same as Figure [T] but clusters are binned in Mvir- 




1 10 

Figure 9. The mass-concentration relation obtained from the 
stacked lensing analysis. We show stacking results of 3 6^ bins 
{filled squares) and 3 Mvir bins (open circles). The mass and con- 
centration measured from stacked strong and weak lensing a nal- 
ysis of 4 massive clusters at z ~ 0.32 ( Umetsu et al.llioilbl ) are 
indicated by a cross. The black shading region shows theoretically 
expected mass-concentration relations with the lensing bias (see 
Appendix |A] for details). The solid and dotted lines are bets-fit 
relation from individual analysis shown in Figure [5] 
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Figure 10. Concentration parameters from stacking lensing anal- 
ysis as a function of the Einstein radius for the source redshift 
Zs = 2. The black shading region indicates the theoretical ex- 
pectation with the selection effect (see Appendix IA] for details). 
Symbols are same as Figure [9] 



fit model predicts (^e = 9^1), although they are consis- 
tent with each other well within la. The best-fit mass and 
concentration are (Mvir) = 4.57to 3i x IO^^H'^Mq and 
(cvir) = 5.75^0 57- note that the mean mass measured 
by the stacking analysis agrees well with the mean mass of 
strong lens selected clusters predic ted by ray-tracing; si mu- 
lati ons, (Mvir) ^ 4.2 X 10 ^^ h~^M(:^ (iHennawi et al.ll2007l . see 
also iBavliss et aLll2Qllbl ) . 

The concentration parameter measured in the stacked 
tangential shear profile is broadly consistent with the re- 
sult of individual analysis (see Figure [S} , but appears to be 
slightly smaller than the mass-concentration relation con- 
strained from our lensing sample. Here we estimate the im- 



pact of the possible averaging effect by computing an average 
shear profile from reduced shear profiles of the NFW profile 
with best-fit values of the mass and concentration from indi- 
vidual cluster analysis (Table [4|), and comparing it with the 
reduced shear profile of the NFW profile with the median 
values of the mass and concentration. We find that both 
profiles agree well near the virial radius, but the averaged 
profile underestimates the shear profile toward the centre, 
with ^ 10% systematic difference at around 0.2/i~^Mpc. 
Both the profiles agree well again at the innermost radii 
of < 0.1/i~^Mpc down to the strong lensing region. The 
systematic difference translates into the concentration pa- 
rameter of Acvir —1.2, and hence it can partly explain 
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the smaller concentration parameter value from the stacked 
lensing analysis. 

Another possible reason for the smaller concentration 
from the stacked tangential shear is a wide range of of our 
sample, which results in the large error on the mean Einstein 
radius and therefore in the much weaker constraints from 
strong lensing compared with individual modelling cases. 
Hence, we conduct the same stacking analysis by dividing 
our cluster sample into 3 bins. In order to assure reason- 
able constraints from strong lensing, we remove 9 clusters 
which have large errors on 0^ mostly because of the lack of 
arc redshift information. To test the mass dependence of the 
concentration, in addition to ^e bins we consider 3 Mvir bins 
too. We use all the 25 clusters for the mass bin analysis. 

Results of our stacking analysis in different ^e and Mvir 
bins are shown in Figures [71 and (8] respectively, and are 
summarized in Table O We find that clusters in the largest 
^E bin are indeed most massive. However, the second and 
third ^E bins have similar mean virial masses, and the dif- 
ference of the Einstein radii appear to be derived by the 
different concentrations. On the other hand, different mass 
bins have similar Einstein radii, but the concentrations are 
clearly larger for smaller masses. 

These results can be used to check the mass- 
concentration relation inferred from individual analysis of 
clusters. Figure [9] shows the mass-concentration relation 
similar to Figure [S] but this time the relation obtained 
from stacked lensing analysis. We find that the mass- 
concentration relation from stacking analysis is in reasonable 
agreement with the best-fit relation constrained from indi- 
vidual analysis of strong and weak lensing (equation [26)) . In 
particular, the strong mass dependence of the concentration 
is clearly seen in the stacking analysis as well. The slightly 
smaller normalization compared with individual analysis can 
partly be ascribed to the averaging effect as described above. 
Thus the stacking analysis further confirms the measure- 
ment of the mass-concentration relation from our sample of 
clusters. 

In Figure [10] we study the dependence of concentration 
parameters derived from the stacking analysis with the Ein- 
stein radius. In particular we compare it with semi-analytic 
calculation conducted in Appendix [A] which predicts that 
the clusters with larger Einstein radii are more concentrated. 
While it is hard to see this trend in our cluster sample, we 
find that the hi gh conce ntr ation o f massive lensing clusters 
presented bv fUmetsu et aP (|2011bl) can be explained in this 
conte xt. Our result suggests that the average concentration 
of the lUmetsu et al.1 (|2011bl ) cluster sample is in good agree- 
ment with the theoretical expectation given the very large 
Einstein radii of 6'e ~ 40^^ . 

5.2 Two-dimensional stacking analysis 

In addition to stacking of the tangential shear profile, we 
conduct stacking of two-dimensional (2D) shear maps to 
study the mean shape of the projected dark matter dis- 
tribution in clusters. Such 2D stacking analysis has been 
attem p ted for samples of ga la xies (Nat araian & Ref regier 



20001: iHoekstra et al.1 l2004l : iMandelbaum et al l2006al : 
Parker et al.1 1200 /t ) or for clusters (jEvans Bridy l2009h . 
The biggest problem of these 2D stacking analysis has been 
that the position angle (orientation) of the projected mass 



distribution has to be known for each cluster when stacking. 
In previous work it was assumed that the position angle of 
the mass distribution coincides with that of the light distri- 
bution, e.g., the surface brightness distribution of the cen- 
tral or satellite galaxy distributions , although the assump- 
tion has not yet been fully justified Oguri et al. l l2010l : lBetd 

Honl). 

Our unique sample of strong and weak lensing clusters 
provides an important means of overcoming this difficulty. 
The idea is that strong lens modelling can generally con- 
strain the position angle of the dark halo component quite 
well, which can be used as a prior information for the po- 
sition angle to stack weak lensing signals. This procedure 
evades any assumptions on the alignment between mass and 
light, and hence should enable much more robust 2D stack- 
ing analysis. 

As in Section 15. 1[ we conduct stacking analysis in the 
physical length scale. For each cluster, we adopt the position 
angle obtained in strong lens mass modelling {Oe in Table [2]) 
to rotate the catalogue of the background galaxies by —Oe 
such that the the position angle of the dark halo is aligned 
with the North-South axis. Specifically, the position of a 
galaxy at (x, y) with respect to the cluster centre is changed 
to 

X — X cos + 2/ sin , (28) 
y' = — xsin^e + 2/cos^e, (29) 

and the two shear components (^i, ^2) are modified as 

g'l — 5'i cos 2^e + 5'2 sin 2^e, (30) 
§2 — —gism20e + g2cos20e. (31) 

We stack the rotated shear catalogue in the physical unit, 
Ticrg'{r'), to obtain the average 2D shear map of our cluster 
catalogue. The cluster catalogues analyzed in this section is 
same as those in Section [5.11 containing 25 clusters in total. 

The stacked 2D shear map, as well as the corresponding 
density map reconstructed from the shear map, are shown 
in Figure [TT1 As expected, the projected mass distribution 
from the stacked 2D shear map is quite elongated along the 
North- South direction, suggesting the highly elongated mass 
distribution of our cluster sample. As a sanity check, we also 
compute the 2D shear map without any alignment of the 
position angle when stacking. The resulting mass distribu- 
tion shown in Figure [TTI appears to be circular symmetric, 
which supports that the highly elongated distribution in our 
stacked map is not an artifact. 

We constrain the ellipticity of the projected 2D mass 
distribution by directly fitting the 2D shear map with the 
elliptical NEW model prediction. Here w e closely follow the 
procedure detailed in I Oguri et~aL for the 2D shear 

fitting. Briefiy, we modify the convergence hi{r) (i.e., the 
projected surface mass density) of the spherical NEW pro- 
file simply by introducing the ellipticity in the iso-density 
contour as ^^/(l ~ c) + 2/^(1 — e). With this proce- 

dure our definition of the ellipticity ise=l — 6/a, where 
a and b are major and minor axis lengths of the isoden- 
sity contour. The corresponding shear pattern is computed 
by solving the Poisson equation. We then construct pix- 
elized distortion field by computing mean shear and errors 
in each bin, and compare it with the elliptical model pre- 
diction, adopting the pixel size of 0.1/i"^Mpc. We add the 
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Figure 11. The two-dimensional weak lensing shear maps obtained from stacking analysis of 25 clusters. The sticks shows observed 
directions and strengths of weak lensing shear distortion. Colou r contours are the surface density map reconstructed from the shear map 
using the standard inversion technique (|Kaiser Squires|[l993D . Both the shear and density maps are smoothed for illustrative purpose. 
Left: The result when the position angle of each cluster is aligned to the North-South axis before stacking, by using the position angle 
measured in strong lens modelling. The resulting stacked density distribution is clearly elongated along the North-South direction. Right: 
The result without any alignment of the position angle when stacking. The resulting density distribution is nearly circular symmetric in 
this case. 



contribution of t he large scale stru cture to the error covari- 
ance matrix fsee lOguri et al.|[2010[ ). We perform fitting in a 
6/i~^Mpc X 6/i~^Mpc region, but remove the innermost 4x4 
pixels considering several possible systemati cs that might b e 
affecting signals near the centre. Unlike Q guri et alJ (|2010l ). 
we fix the mass centre to the assumed centre (the position 
of the brightest galaxy in strong lensing region), because 
strong lensing available for our cluster sample allows a reli- 
able identification of the mass centre for each cluster. Thus 
we fit the 2D shear map with four parameters (Mvir, Cvir, e, 
^e), employing a Markov Chain Monte Carlo technique. 

In Figure [T2I we show the posterior likelihood distribu- 
tion of the mean ellipticity (e) from the 2D stacking analysis 
of all the 25 clusters. When the position angles are aligned, 
the resulting density distribution is indeed quite elliptical 
with the mean ellipticity of (e) = 0.47 =b 0.06. We find that 
the elliptical NFW model improve fitting by Ax^ = 26.9 
compared with the case e = 0, which indicates that the 
detection of the elliptical mass distribution is significant 
at the 5(7 level. The measured mean ellipticity is consis- 
tent with the average ellipticity from strong lens modelling 
(e) = 0.38 =b 0.24, although the latter involves large scat- 
ter. The best-fit position angle of Oe — 9- 1^4.1 deg slightly 
deviates from the expected position angle of = 0, but 
they are consistent with each other within 2a (Ax^ < 4). In 
contrast, if the position angles are not aligned in stacking 
shear signals, the resulting constraint on the mean elliptic- 
ity is (e) < 0.19, i.e., it is fully consistent with the circular 
symmetric mass distribution e — within la. 

We compare this result with the theoretical prediction 
in the A CDM model. For th is purpose we employ a triaxial 
model of Ijing Siltol (I2OO2I I. Assuming that the halo orien- 



Table 6. Summary of the two-dimensional stacking analysis 
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^E-3 


29+°-^^ 
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tation is random, we compute the probability distribution 
of the ellipticity by proj ecting the triaxia l halo along arb i- 
trary directions (jOguri e t al. 2003; Qguri KeetonI l2004l l . 
In this analysis we fix the mass and redshift of the halo to 
Mvir = 4.6 X lO^'^/i~^M0 and z = 0.469, which are mean 
mass and redshift of the 25 clusters. We find that the mean 
ellipticity predicted by this model is (e) = 0.44, in excellent 
agreement with the measured ellipticity. The analysis pre- 
sented in Appendix [A] indicates that the effect of the lensing 
bias on the mean ellipticity is small, with a possible shift of 
the mean ellipticity of < 0.05 at most, and hence it does not 
affect our conclusion. Our result is also in good agreement 
with the pre vious lensing measurement of the ellipticity by 
lOguri et al.1 (|201Q| ) in which 2D shear maps of individual 
clusters are fitted with the elliptical NFW profile, rather 
than examining the stacked shear map. 

We check the sensitivity of our ellipticity result on the 
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Figure 12. The marginalized probability distribution of the 
mean ellipticity (e) from stacked weak lensing analysis of 25 
clusters. The solid line indicates the case when the position an- 
gles are aligned according to those measured with strong lens 
modelling, in which the mean ellipticity is detected at 5cr level 
((e) = 0.47 ±0.06). The dashed line is the marginalized probabil- 
ity distribution for stacking without any alignment of the position 
angles, for which the mass distribution is consistent with the cir- 
cular symmetric distribution ((e) < 0.19). The vertical dotted 
line indicates the theo retical expectation, (e) = 0.44, based on a 
triaxial halo model of Ijing Sutol (|2002l ). The open circle with 
errorbar shows the average ellipticity and la scatter from strong 
lens modelling. 



size of the fitting region, as one possible concern is that 
infalling matter associated with the filamentary structure 
outside clusters might boost the mean ellipticity. Figure [13] 
shows how the constraint on the mean ellipticity changes 
by making the size of the fitting region smaller from our 
fiducial choice (half the box size of 3/i"^Mpc). The Fi gure 
indicates that the detection of the mean ellipticity of ^ 0.45 
is robust against the choice of the fitting size, as the results 
are consistent down to half the box size of ^ 0.3/i~^Mpc 
where the constraint become significantly weak. The anal- 
ysis also implies that the ellipticity of the projected mass 
distribution does not change very much with radius. Theo- 
retically, da rk haloes are expe cted to be more elongated near 
the centre ( Jing Sutdl20o3 l , although the effect of baryon 
cooling and star formatio n can m ake the shape rounder par- 
ticula rly near the centre (jKazantzidis et al.ll2004l : iLau et al.l 
l201lh . Our detection of the elliptical mass distribution down 
to small radii may therefore help constraining the amount 
of cooling in clusters. 

Finally we check the mean ellipticities for different 
and Mvir bins. Figure [T^ shows the results, and Table [6] gives 
the summary of constraints on the mean ellipticities. We find 
a possible decrease of (e) at the largest Einstein radius bin, 
although it is not very significant given the large errorbars. 
On the other hand, there appears no simple trend in the 
mass bin result. 

While results for all subsamples are almost con- 
sistent with the full sample result, the dependence on 

might be suggestive of conflicting selection effects 
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Figure 13. The mean ellipticities measured in stacked 2D shear 
map as a function of the box size for fitting. Our fiducial result and 
its 1(7 range adopting half the box size of 3/i~-'^Mpc are indicated 
by the horizontal line with shading. 



on our sample. There are two main lensing biases that 
can affect the results on the mean ellipticity measure- 
ment. One is the orientation bias, i.e., clusters with 
larger Einstein radii appear to be rounder than nor- 
mal clusters because of the alignm ent of the major 
axis with the line-of-sight direction (|Oguri BlandfordI 
The other selection effect is directly related to 
the ellipticity of the projected mass distribution. Because 
the ellipticity significantly enhances the giant arc cross 
section (e.g., .Meneghetti/Bartelmann. Moscard"iiiil l2003l : 



iMeneghetti et al.l 2007 ). a sample of clusters with promi- 
nent lensed arcs should have more elliptical projected mass 
distribution. These two selection effects apparently conflict 
with each other, although a naive expectation is that the 
orientation bias dominates for large Einstein radii, whereas 
the ellipticity bias is more significant for small Einstein 
radii. Indeed, our semi-analytic calculations presented in 
Appendix [Al shows a clear dependence of the mean ellipticity 
on the Einstein radius, as plotted in Figure [T^ The observed 
trend appears to be consistent, at least qualitatively, with 
the theoretical prediction. A large sample is needed to con- 
firm this trend more robustly. 



6 RELATION BETWEEN MASS AND 
MEMBER GALAXY DISTRIBUTIONS 

Understanding the relation between the mass and clus- 
ter member galaxy distributions is importa nt for the cos- 
mological use of optica l clusters (e-g - ; iKoester et aDl2007l : 
IRozo et all l2009l . l201ll : iRvkoff et all l201lh as well as de- 
t ailed stacking analy sis for a large sample of clusters (e.g., 
lEvans k B^idl3 12009 ^1 . The spatial distribution of galaxies, 
in comparison with that of subhaloes in an A/'-body simula- 
tion, is crucial for explorin g the formation history o f clus- 
ter member galaxies fe.g., lNagai Kravtsovl l2005l V Here 
we examine spatial distributions of cluster member galax- 
ies using the stacking technique, and compare them with 
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Table 7. analysis results on the member galaxy distribution 
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Figure 14. Mean ellipticities from stacked shear maps in three 
different (upper panel) and Mvir (lower panel) bins. The hor- 
izontal line with shading indicate best-fit and la range of the 
full sample result. Dashed lines show the predicted mean ellip- 
ticity from semi-analytic calculations with the lensing bias (see 
Appendix |A] for details). 



accurate measurements of mass distributions from stacked 
lensing analysis. 

While the colour information from ^ri-band images al- 
lows us to select the galaxy population around a given red- 
shift efficiently, the selection is not perfect in the sense 
that there is some contamination from foreground and back- 
ground galaxies. Thus, in this paper we take advantage of 
the high number density of galaxies after stacking to sub- 
tract the foreground and background galaxies statistically. 
Specifically, we assume that the stacked number density dis- 
tribution is described by the sum of cluster member galaxies 
and background galaxies as: 



(r) + i 



(32) 

Snfw is the surface density 



where ag is the normalization, 
distribution of the NFW profile, and bg is the number den- 
sity of foreground/background galaxies which is assumed to 
be constant across the field. 

In this paper, we focus on two parameters that describe 
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the number density distribution of member galaxies. One is 
the scale radius rs of the NFW profile, which is used to check 
whether or not the galaxy distribution is more extended than 
the mass distribution. The other is the ellipticity of the pro- 
jected number density distribution, as our stacking analysis 
enables detailed analysis of 2D distributions for both mass 
and member galaxy distributions. We constrain the scale 
radius rs by fitting azimuthally averaged radial profile of 
the galaxy number density, whereas the ellipticity is derived 
from the 2D fitting of the density distribution. The total 
number of parameters for the radial profile fitting is 3 (rs, 



ttp, and bg), and the that of 2D fitting is 5 (e, Oe, Ts, a^, and 

We consider two distinct galaxy populations, one is red 
member galaxies selected by the colour cut (equations [71 and 
(8]), and the other is all galaxies without any colour cut. For 
each galaxy population, we consider galaxies that are more 
luminous than the luminosity cut, i.e., L > Lent- The lu- 
minosity cut is defined in terms of L*, where L* is the lu- 
minosity corresponding to the z-band abso lute magnitude 
Mt = -21.22 + 51og/i (jRvkoff et al.ll201lh . For the high- 
est cluster redshift in our sample, a luminosity O.IL* corre- 
sponds to i ^ 24.5, which is well above the magnitude limit 
of our Subaru imaging. 

We show our fitting results in Figure [15] and Table (T] We 
find that the mass and galaxy number distributions agree 
reasonably well for the luminous galaxies, L > L*, for both 
the scale radius and ellipticity. However, for fainter galaxies 
of L ^ O.IL* the scale radius of the galaxy distribution tends 
to be larger (i.e., smaller concentration parameter values), 
and the average ellipticity is smaller. 

The radial distribution has also been studied in pre- 
vious work. Measurements of the projected number den- 
sity distributions suggest low concent ration values o f c ~ 
(Lin, Mohr Stanford 2004: 



Hansen et al.ll2005l: iLin k Mohrll200^ : 

3 



Katgert et al.1 12004 



Biviano &: Poggiantil 



20091 ) . iLin k Mohri ^200'^ found that more luminous galax- 



ies tend to be distributed with large concentrations, which is 
consistent with our results in Figure [T5l Simulation work by 
.Nagai Kravtsovl (j2005) also predicts that the galaxy num- 
ber density distribution is slightly more extended than the 
mass distribution, which is again consistent with our result, 
although no significant dependence of the concentration on 
the stellar mass was reported from the simulations. 

There have been measurements of the ellipticity of 
the projected galaxy number density distri butions (e.g., 
IPaz et al.|[2006l : iNie derste-O stholt et al.l l2010V but its con- 
nection with the underlying mass distribution has not yet 
fully been explored. The similar ellipticity for luminous 
member galaxies is reasonable, but the origin of the smaller 
mean ellipticity for fainter galaxies is unclear. Essentially, 
the smaller mean ellipticity for the stacked cluster sample 
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Figure 15. The mean scale radius (r^) {upper) and mean elliptic- 
ity (e) {lower) of the member galaxy distribution from stacking 25 
clusters, as a function of the luminosity cut Lcut- Filled squares 
and open circles show the results for red member galaxies and all 
member galaxies, respectively. Horizontal lines with shading are 
the scale radius and ellipticity for the mass distribution, measured 
by stacked lensing, for the same cluster sample. 



indicates either the galaxy distribution is on average rounder 
or there is a large misalignment (i.e., difference in the po- 
sition angles) between the mass and galaxy distribution. In 
either case careful simulation work is crucial to check if the 
observed distribution is in agreement with current under- 
standing of cluster galaxy formation. 



7 CONCLUSION 

We have performed a combined strong and weak lensing 
analysis for a sample of 28 clusters in the redshift range 
0.28 < z < 0.68. The cluster sample is based on the SGAS, a 
large survey of giant arcs amongst SDSS clusters. In this pa- 
per, we have presented extensive follow-up observations with 
Subaru/Suprime-cam, which enables reliable weak lensing 
measurements out to large radii. Combined with the giant 
arcs with significant amount of spectroscopic information, 
we can constrain mass distributions of these clusters quite 



well to measure the concentration and shape of the mass 
distribution. 

The mass-concentration relation derived from our lens- 
ing analysis has the slope of Cvir cx: M^^"^^^°"^^, which is 
significantly steeper than the slope predicted by the theoret- 
ical expectation of —0.2 which includes the effect of the lens- 
ing bias. Concentrations measured by the combined lensing 
analysis are in reasonable agreement with the ACDM pre- 
diction for massive clusters (Mvir lO^^/i-^M©), if we take 
proper account of the lensing bias. The result indicates that 
the anoma lously h i gh concentration jB roadhurst et al.'' 2005l, 
2008; Qgu ri et al.1 l2009l : lUmetsu et a l. 2011b; Gralla eTahl 
[20111), which has been claimed from analysis of small num- 
ber of lensing clusters, is now much less evident, thanks to 
the much larger number of clusters we have analyzed in the 
paper. On the other hand, observed concentrations appear 
to be significantly higher than theoretical expectations for 
lower masses of Mvir - lO^^/i-^Mo, which may be explained 
as arising from the effects of cooling baryons in cluster cen- 
tres. 

We have also stacked tangential shear profiles for the 
sample of clusters to obtain an accurate mean profile, which 
is seen to be in good agreement with the NFW profile. Our 
stacking analysis at different Einstein radii and mass bins 
has confirmed our results on the mass-concentration rela- 
tion from the individual analysis of these clusters. In addi- 
tion to the radial profile, we have explored the stacked 2D 
shear map to study the shape of the projected mass distribu- 
tion. By aligning shear maps of individual clusters with the 
position angles of dark haloes from strong lens modelling, 
we were able to detect the elliptical shape of the stacked 
mass distribution at the 5cr level. The mean ellipticity of 
(e) = 0.47 ± 0.06 is in excellent agreement with the ACDM 
expectation. The significant detection of the highly elliptical 
shape of dark matter hal oes wi th weak lensing has also been 
reported bv lOguri et al.l (|201Qh : this work confirms the pre- 
vious finding using a different sample and technique. Finally, 
based on stacking analysis we have compared distributions 
of cluster member galaxies with mass distributions, finding 
good agreement between them for luminous (L > L*) mem- 
ber galaxies. Distributions of fainter galaxies are found to 
be more extended. 

Our work has demonstrated the power of combined 
strong and weak lensing, not only for analysis of individ- 
ual clusters but also for detailed stacked lensing analysis. 
The techniques described in the paper can be applied to 
the unique sample of lensing clusters obtained by the Clus- 
ter Lensing And Su pernova survey with Hubble (CLASH; 
Postman et al. I l201lh . The detailed comparison of SGAS re- 
sults as presented in the paper with upcoming CLASH re- 
sults, as well as lensing resu lts from the Local Cluster Sub- 
structure Survey fLoCuSS: iOkabe et al.ll2010l ). will be very 
useful, particularly because of markedly different sample se- 
lection between these surveys. 
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APPENDIX A: EXPECTED PROPERTIES OF 
STRONG LENSING SELECTED CLUSTERS 

We predict properties of clusters in our strong-lens selected 
cluster sam ple using a semi-analytic model developed by 
lOguri &; Bl andford (2009), which is based on a triaxial halo 
model of [jing Sutd ((2003) • In brief, a catalogue of haloes 
are generated according to the mass function and axis ra- 
tio distribution derived from A^-body simulations, and each 
halo is projecte d along ra ndom direction to compute its lens- 
ing property (O guri et al. 2003; Oguri & Keeton 2004). The 
projected convergence profile is compared with that of a 
spherical NFW profile to estimate the mass Mvir,2D, the 
concentration Cvir,2D, and the ellipticity e of the projected 
mass distribution. Here we fix the lens (cluster) and source 
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redshift to zi — 0.45 and z — 2^ respectively, which are typ- 
ical for our cluster lens sample analyzed in the paper. 

We compute a giant arc cross section (Jarc of each halo 
from ray-tracing of extended sources using glafic (Q guril 
The source is assumed to have a Sersic profile with 
Te = 0^.^4 and the ellipticity randomly assigned between 
and 0.5. For each cluster we compute the arc cross section 
for the length-to- width ratio oil/w > 5 by randomly throw- 
ing source galaxies in the source plane. Output images are 
convolved with the Gaussian kernel with the FWHM of 0^8 
in order to take account of the seeing effect. 

We study the impact of the lensing bias by averaging 
concentrations and ellipticities of the halo catalogue with an 
appropriate weight that mimics our selection criteria. Obvi- 
ously our cluster sample is weighted by the arc cross section 
(Tare, which cau be a reasonable choice of the weight. Fur- 
thermore, we preferentially conduct follow-up observations 
for clusters with larger Einstein radii (i.e., giant arcs located 
more distant from the cluster centre), because they are ex- 
pected to be more massive. We can model this selection ef- 
fect, e.g., by multiplying to the weight. 

Figure lAll shows the mass-concentration relation de- 
rived from the projected mass distribution assuming a spher- 
ical mass distribution, which is relevant for the comparison 
with our lensing measurement, based on the semi- analytic 
calculation described above. We find that the enhancement 
of the concentration parameter due to the lensing bias is a 
strong function of the mass. For instance, when both the 
arc cross section and the Einstein radius are used for com- 
puting the lensing bias, the concentration parameter is en- 
hanced by -80% for Mvir - 8 x 10^^ H'^Mq, -30% for 
Mvir - 4 X lO^^/i-^Mo, and -20% for Mvir - lO^^h'^MQ. 

We also compare our result based on the triaxial model 
with other work studying the mass-concentration relation 
in A^-body simulations based on the spherical NEW pro- 
file. We find that our mass-concentration relatio n has larger 
concentration than predicted by the relation of iDuffv et al.l 
(|2008[ ). which has b een adopted for co mparisons with lens- 
ing measurements i n lOguri et al.1 (|2009'l. On the other hand, 
IPrada et al.1 (|201lh concluded that concentrations of mas- 
sive haloes evolve little with the redshift. As a result, con- 



centrations predicted by the mod el of IPrada et al. (|201lh 
are much larger compared with the iDuffv et aL I (|2008l ) rela- 
tion predicts. We note th at such littl e evolution of massive 
haloes have been noted bv Zhao et al.l (|2003 ) and lZhao et al.l 
(|2009h . Our mass-concentration relation based on the tri- 
axial model resides in between the iDuffv et all (|2008h and 
IPrada et al.1 (l201lh relations, and more resemble the rela- 
tion presented by iGao et aU (|2008 ^1 which was essentially 
the m odification of the model proposed by iNavarro et al.l 
(|l997f ). 

We also investigate the impact of the lensing bias 
on the ellipticity of the projected mass distribution. We 
find that the lensing bias due to the arc cross section 
tends to increase the mean ellipticity, which is under- 
stood by the fact that the ellipticity of the projected 
mass distribution significantl y enhances the giant arc cross 
section (e.g., Men eghetti, Bartelmann, Moscardinil |2003l: 
iMeneghetti et al I I2OO7I I . In contrast, the lensing bias due 
to the Einstein radius decreases the mean Einstein radius, 
because of the alignment of the major ax is with the line-of- 
sight direction (jOguri k Blandfordllioool ). Therefore, these 
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Figure Al. Upper: The mass-concentration relation of the pro- 
jected mass d istribution fo r z = 0.45, predicted by a triaxial 
halo model of Ijing Sutd (|20Q2l V Crosses indicate the average 
concentration parameter from our semi-analytic calculation with- 
out any weighting, whereas filled triangles and open circles are 
average concentration parameters with the lensing bias, assum- 
ing weights from giant arc cross sections (Tare or arc cross sec- 
tions plus the Einstein radii ^e, respectively Lines are predicted 
mass-concentration relation f rom spheric al averaging of halo es in 
various A^-b ody simulations (|Gao et al .> 2008; Duffv et al. 20^; 
IPrada et al.ll201lD . Lower: Similar to the upper panel, but aver- 
age ellipticities of the project mass distribution as a function of 
the mass are shown. 



two lensing biases counteract with each other. In either case, 
however, the effect of the lensing bias appears to be small, 
wit h the change of the m ean ellipticity by ^ 0.05 at most. 

lOguri BlandfordI (^2009. ) has shown that extreme lens- 
ing clusters having very large Einstein radii are more 
severely affected by the lensing bias. Thus we check the de- 
pendence of the concentration and projected ellipticity on 
the Einstein radius ^e- For the mass distribution of the clus- 
ter, we assume a flat prior of log Mvir between 10^^h~^ Mq 
and lO^^/i~^M0, and the number distribution predicted by 
the mass function of dark haloes above 10^^h~^ Mq, which 
is expected to more or less resemble the selection function of 
our cluster sample. The result shown in Figure lA2l indicate 
that both the concentration and ellipticity depend strongly 
on the Einstein radius such that clusters with larger Ein- 
stein radii are more concentrated and spherical. Again, this 
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Figure A2. The average concentration {filled squares) and el- 
lipticity {open triangles) of the projected mass distributions as a 
function of the Ein stein radius ^e- Th ey are derived from a tri- 
axial halo model of IJ ing Sutol (120021) and include the effect of 
the lensing bias from the arc cross section. 



can be interpreted by the alignment of the major a xis with 
the line-of-sight direction (|Oguri k, BlandfordlbOOOh . 



APPENDIX B: LENSING ANALYSIS FOR 
INDIVIDUAL CLUSTERS 

We show for each cluster the critical curve of the best-fit 
strong lens modelling plotted on the Subaru/Suprime-cam 
^fri-composite image {2' x 2'] squares are positions of multi- 
ple images used for mass modelling) , weak lensing mass map 
plotted on the Subaru/Suprime-cam r-band image (contours 
are drawn with spacing of la noise level, and the cross indi- 
cates the position of the brightest cluster galaxy), and the 
tangential shear profile as well as the best-fit NFW profiles. 
The shading in the tangential profile plot indicates the mea- 
sured Einstein radius and its la error at the arc redshift. 
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Figure Bl. SDSSJ0851+3331, SDSSJ0915+3826, SDSSJ0957+0509, SDSSJ1004+4112 
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Figure B2. SDSSJ1029+2623, SDSSJ1038+4849, SDSSJ1050+0017, RCS2J1055+5548 
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Figure B3. SDSSJlllO+6459, SDSSJ1115+5319, SDSSJ1138+2754, SDSSJ1152+3313 
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Figure B4. SDSSJ1152+0930, SDSSJ1209+2640, SDSSJ1226+2149, SDSSJ1226+2152 
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Figure B5. A1703, SDSSJ1315+5439, GHO132029+3155, SDSSJ1329+2243 
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Figure B6. SDSSJ1343+4155, SDSSJ1420+3955, SDSSJ1446+3032, SDSSJ1456+5702 
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Figure B7. SDSSJ1531+3414, SDSSJ1621+0607, SDSSJ1632+3500, SDSSJ2111-0114 
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